Exploration of cell type-specific somatic mutations in schizophrenia and the impact of maternal immune activation on the somatic mutation profile in the brain.

AIM: Schizophrenia (SZ) is a severe psychiatric disorder caused by the interaction of genetic and environmental factors. Although somatic mutations that occur in the brain after fertilization may play an important role in the cause of SZ, their frequencies and patterns in the brains of patients and related animal models have not been well studied. This study aimed to find somatic mutations related to the pathophysiology of SZ. METHODS: We performed whole-exome sequencing (WES) of neuronal and nonneuronal nuclei isolated from the postmortem prefrontal cortex of patients with SZ (n = 10) and controls (n = 10). After detecting somatic mutations, we explored the similarities and differences in shared common mutations between two cell types and cell type-specific mutations. We also performed WES of prefrontal cortex samples from an animal model of SZ based on maternal immune activation (MIA) and explored the possible impact of MIA on the patterns of somatic mutations. RESULTS: We did not find quantitative differences in somatic mutations but found higher variant allele fractions of neuron-specific mutations in patients with SZ. In the mouse model, we found a larger variation in the number of somatic mutations in the offspring of MIA mice, with the occurrence of somatic mutations in neurodevelopment-related genes. CONCLUSION: Somatic mutations occurring at an earlier stage of brain cell differentiation toward neurons may be important for the cause of SZ. MIA may affect somatic mutation profiles in the brain.

LINE-1 hypomethylation, increased retrotransposition and tumor-specific insertion in upper gastrointestinal cancer.

The long interspersed nuclear element-1 (LINE-1) retrotransposons are a major family of mobile genetic elements, comprising approximately 17% of the human genome. The methylation state of LINE-1 is often used as an indicator of global DNA methylation levels and it regulates the retrotransposition and somatic insertion of the genetic element. We have previously reported the significant relationship between LINE-1 hypomethylation and poor prognosis in upper gastrointestinal (GI) cancers. However, the causal relationships between LINE-1 hypomethylation, retrotransposition, and tumor-specific insertion in upper GI cancers remain unknown. We used bisulfite-pyrosequencing and quantitative real-time PCR to verify LINE-1 methylation and copy number in tissue samples of 101 patients with esophageal and 103 patients with gastric cancer. Furthermore, we analyzed the LINE-1 retrotransposition profile with an originally developed L1Hs-seq. In tumor samples, LINE-1 methylation levels were significantly lower than non-tumor controls, while LINE-1 copy numbers were markedly increased. As such, there was a significant inverse correlation between the LINE-1 methylation level and copy number in tumor tissues, with lower LINE-1 methylation levels corresponding to higher LINE-1 copy numbers. Of particular importance is that somatic LINE-1 insertions were more numerous in tumor than normal tissues. Furthermore, we observed that LINE-1 was inserted evenly across all chromosomes, and most often within genomic regions associated with tumor-suppressive genes. LINE-1 hypomethylation in upper GI cancers is related to increased LINE-1 retrotransposition and tumor-specific insertion events, which may collectively contribute to the acquisition of aggressive tumor features through the inactivation of tumor-suppressive genes.

Quantification of cytosine modifications in the aged mouse brain.

Quantifying cytosine modifications in various brain regions provides important insights into the gene expression regulation and pathophysiology of neuropsychiatric disorders. In this study, we quantified 5-methylcytosine (5-mC), 5-hydroxymethylation (5-hmC), and 5-formylcytosine (5-fC) levels in five brain regions (the frontal lobe, cerebral cortical region without frontal lobe, hippocampus, basal ganglia, and the cerebellum) and the heart at three developmental periods (12, 48, and 101 weeks). We observed significant regional variations in cytosine modification. Notably, regional variations were generally maintained throughout development, suggesting that epigenetic regulation is unique to each brain region and remains relatively stable with age. The 5-mC and 5-hmC levels were positively correlated, although the extent of the correlations seemed to differ in different brain regions. On the contrary, 5-fC levels did not correlate with 5-mC or 5-hmC levels. Additionally, we observed an age-dependent decrease in 5-fC levels in the basal ganglia, suggesting a unique epigenetic regulation mechanism. Further high-resolution studies using animal models of neuropsychiatric disorders as well as postmortem brain evaluation are warranted.

Sex-dependent behavioral alterations in a poly(I:C)-induced maternal immune activation mouse model without segment filamentous bacteria.

Maternal immune activation is one of the environmental risk factors for offspring to develop psychiatric disorders. A synthetic viral mimetic immunogen, polyinosinic-polycytidylic acid (poly(I:C)), is used to induce maternal immune activation in animal models of psychiatric disorders. In the mouse poly(I:C) model, the existence of segment filamentous bacteria (SFB) in the maternal intestine has been reported to be important for the induction of ASD-related behavioral alterations as well as atypical cortical development called cortical patches. This study aimed to elucidate the effect of a single poly(I:C) injection during embryonic day (E) 9 to E16 on offspring's behavior in the ensured absence of maternal SFB by vancomycin drinking in C57BL/6N mice. The cortical patches were not found at either injection timings with poly(I:C) or PBS vehicle, tested in male or female offspring at postnatal day 0 or 1. Prepulse inhibition was decreased in male adult offspring most strongly at poly(I:C) injection timings later than E11, whereas a modest but significant decrease was observed in female offspring with an injection during E12 to E15. The decrease in social interaction was observed in female offspring most conspicuously at injection timings later than E11, whereas a significant decrease was observed in male offspring with an injection during E12 to E15. In conclusion, this study indicated that behavioral alterations could be induced without maternal SFB. The effect on behavior was substantially different between males and females.

Identification of epigenetically active L1 promoters in the human brain and their relationship with psychiatric disorders.

Long interspersed nuclear element-1 (LINE-1, L1) affects the transcriptome landscape in multiple ways. Promoter activity within its 5'UTR plays a critical role in regulating diverse L1 activities. However, the epigenetic status of L1 promoters in adult brain cells and their relationship with psychiatric disorders remain poorly understood. Here, we examined DNA methylation and hydroxymethylation of the full-length L1s in neurons and nonneurons and identified "epigenetically active" L1s. Notably, some of epigenetically active L1s were retrotransposition competent, which even had chimeric transcripts from the antisense promoters at their 5'UTRs. We also identified differentially methylated L1s in the prefrontal cortices of patients with psychiatric disorders. In nonneurons of bipolar disorder patients, one L1 was significantly hypomethylated and showed an inverse correlation with the expression level of the overlapping gene NREP. Finally, we observed that altered DNA methylation levels of L1 in patients with psychiatric disorders were not affected by the surrounding genomic regions but originated from the L1 sequences. These results suggested that altered epigenetic regulation of the L1 5'UTR in the brain was involved in the pathophysiology of psychiatric disorders.

A Method for Detection of Somatic LINE-1 Insertions at the Single-Cell Level from Postmortem Human Brain.

Long Interspersed Element-1 (LINE-1, L1) is a retrotransposon that has the ability to amplify its copy in the genome autonomously. L1Hs is a human-specific active subtype of L1 reported to amplify its copy in neural progenitor cells causing genomic mosaicism. This chapter describes a new method named NECO-seq (Novel Elements Concentrated-sequence) to identify the genomic locus of L1Hs insertions at the single-cell level. This protocol contains the steps of (1) preparation of neuronal cell nuclei from a postmortem human brain, (2) whole genome amplification from single neural nuclei (snWGA), (3) single nucleotide polymorphisms (SNPs) genotyping for quality control of snWGA products, (4) library preparation for next-generation sequencing to enrich the genomic locus of L1Hs insertions, and (5) bioinformatic analysis to detect novel somatic L1Hs insertions. This method can detect approximately 97% of L1Hs originally existing in reference human genome and approximately 10-20 newly inserted L1Hs copies in a neuronal cell of a postmortem human brain.

Cell-type-specific DNA methylation analysis of the frontal cortices of mutant Polg1 transgenic mice with neuronal accumulation of deleted mitochondrial DNA.

Bipolar disorder (BD) is a severe psychiatric disorder characterized by repeated conflicting manic and depressive states. In addition to genetic factors, complex gene-environment interactions, which alter the epigenetic status in the brain, contribute to the etiology and pathophysiology of BD. Here, we performed a promoter-wide DNA methylation analysis of neurons and nonneurons derived from the frontal cortices of mutant Polg1 transgenic (n = 6) and wild-type mice (n = 6). The mutant mice expressed a proofreading-deficient mitochondrial DNA (mtDNA) polymerase under the neuron-specific CamK2a promoter and showed BD-like behavioral abnormalities, such as activity changes and altered circadian rhythms. We identified a total of 469 differentially methylated regions (DMRs), consisting of 267 neuronal and 202 nonneuronal DMRs. Gene ontology analysis of DMR-associated genes showed that cell cycle-, cell division-, and inhibition of peptide activity-related genes were enriched in neurons, whereas synapse- and GABA-related genes were enriched in nonneurons. Among the DMR-associated genes, Trim2 and Lrpprc showed an inverse relationship between DNA methylation and gene expression status. In addition, we observed that mutant Polg1 transgenic mice shared several features of DNA methylation changes in postmortem brains of patients with BD, such as dominant hypomethylation changes in neurons, which include hypomethylation of the molecular motor gene and altered DNA methylation of synapse-related genes in nonneurons. Taken together, the DMRs identified in this study will contribute to understanding the pathophysiology of BD from an epigenetic perspective.

Comprehensive DNA Methylation Analysis of Human Neuroblastoma Cells Treated With Haloperidol and Risperidone.

Accumulating evidence suggests that the epigenetic alterations induced by antipsychotics contribute to the therapeutic efficacy. However, global and site-specific epigenetic changes by antipsychotics and those shared by different classes of antipsychotics remain poorly understood. We conducted a comprehensive DNA methylation analysis of human neuroblastoma cells cultured with antipsychotics. The cells were cultured with low and high concentrations of haloperidol or risperidone for 8 days. DNA methylation assay was performed with the Illumina HumanMethylation450 BeadChip. We found that both haloperidol and risperidone tended to cause hypermethylation changes and showed similar DNA methylation changes closely related to neuronal functions. A total of 294 differentially methylated probes (DMPs), including 197 hypermethylated and 97 hypomethylated DMPs, were identified with both haloperidol and risperidone treatment. Gene ontology analysis of the hypermethylated probe-associated genes showed enrichment of genes related to the regulation of neurotransmitter receptor activity and lipoprotein lipase activity. Pathway analysis identified that among the DMP-associated genes, SHANK1 and SHANK2 were the major genes in the neuropsychiatric disorder-related pathways. Our data would be valuable for understanding the mechanisms of action of antipsychotics from an epigenetic viewpoint.

Decreased DNA methylation at promoters and gene-specific neuronal hypermethylation in the prefrontal cortex of patients with bipolar disorder.

Bipolar disorder (BD) is a severe mental disorder characterized by repeated mood swings. Although genetic factors are collectively associated with the etiology of BD, the underlying molecular mechanisms, particularly how environmental factors affect the brain, remain largely unknown. We performed promoter-wide DNA methylation analysis of neuronal and nonneuronal nuclei in the prefrontal cortex of patients with BD (N = 34) and controls (N = 35). We found decreased DNA methylation at promoters in both cell types in the BD patients. Gene Ontology (GO) analysis of differentially methylated region (DMR)-associated genes revealed enrichment of molecular motor-related genes in neurons, chemokines in both cell types, and ion channel- and transporter-related genes in nonneurons. Detailed GO analysis further revealed that growth cone- and dendrite-related genes, including NTRK2 and GRIN1, were hypermethylated in neurons of BD patients. To assess the effect of medication, neuroblastoma cells were cultured under therapeutic concentrations of three mood stabilizers. We observed that up to 37.9% of DMRs detected in BD overlapped with mood stabilizer-induced DMRs. Interestingly, mood stabilizer-induced DMRs showed the opposite direction of changes in DMRs, suggesting the therapeutic effects of mood stabilizers. Among the DMRs, 12 overlapped with loci identified in a genome-wide association study (GWAS) of BD. We also found significant enrichment of neuronal DMRs in the loci reported in another GWAS of BD. Finally, we performed qPCR of DNA methylation-related genes and found that DNMT3B was overexpressed in BD. The cell-type-specific DMRs identified in this study will be useful for understanding the pathophysiology of BD.

Identification and functional characterization of the extremely long allele of the serotonin transporter-linked polymorphic region.

SLC6A4, which encodes the serotonin transporter, has a functional polymorphism called the serotonin transporter-linked polymorphic region (5-HTTLPR). The 5-HTTLPR consists of short (S) and long (L) alleles, each of which has 14 or 16 tandem repeats. In addition, the extralong (XL) and other rare alleles have been reported in 5-HTTLPR. Although they are more frequent in Asian and African than in other populations, the extent of variations and allele frequencies (AFs) were not addressed in a large population. Here, we report the AFs of the rare alleles in a large number of Japanese subjects (N = 2894) consisting of two cohorts. The first cohort (case-control study set, CCSS) consisted of 1366 subjects, including 485 controls and 881 patients with psychosis (bipolar disorder or schizophrenia). The second cohort (the Arao cohort study set, ACSS) consisted of 1528 elderly subjects. During genotyping, we identified 11 novel 5-HTTLPR alleles, including 3 XL alleles. One novel allele had the longest subunit ever reported, consisting of 28 tandem repeats. We named this XL28-A. An in vitro luciferase assay revealed that XL28-A has no transcriptional activity. XL28-A was found in two unrelated patients with bipolar disorder in the CCSS and one healthy subject in the ACSS who did not show depressive symptoms or a decline in cognitive function. Therefore, it is unlikely that XL28-A is associated with psychiatric disorders, despite its apparent functional deficit. Our results suggest that unraveling the complex genetic variations of 5-HTTLPR will be important for further understanding its role in psychiatric disorders.

Establishment of Quantitative PCR Assays for Active Long Interspersed Nuclear Element-1 Subfamilies in Mice and Applications to the Analysis of Aging-Associated Retrotransposition.

The retrotransposon long interspersed nuclear element-1 (LINE-1) can autonomously increase its copy number within a host genome through the retrotransposition process. LINE-1 is active in the germline and in neural progenitor cells, and its somatic retrotransposition activity has a broad impact on neural development and susceptibility to neuropsychiatric disorders. The method to quantify the genomic copy number of LINE-1 would be important in unraveling the role of retrotransposition, especially in the brain. However, because of the species-specific evolution of LINE-1 sequences, methods for quantifying the copy number should be independently developed. Here, we developed a quantitative PCR (qPCR) assay to measure the copy number of active LINE-1 subfamilies in mice. Using the assay, we investigated aging-associated alterations of LINE-1 copy number in several brain regions in wild-type mice and Polg+/D257A mice as a model for accelerated aging. We found that aged Polg+/D257A mice showed higher levels of the type GfII LINE-1 in the basal ganglia than the wild-type mice did, highlighting the importance of assays that focus on an individual active LINE-1 subfamily.

Use of the Illumina EPIC methylation array for epigenomic research in the crab-eating macaque (Macaca fascicularis).

BACKGROUND: Commercially available Illumina DNA methylation arrays (HumanMethylation 27K, HumanMethylation450, and MethylationEPIC BeadChip) can be used for comprehensive DNA methylation analyses of not only the human genome but also other mammalian genomes, ranging from those of nonhuman primates to those of rodents. However, practical application of the EPIC array to the crab-eating macaque has not been reported. METHODS: Through bioinformatic analyses involving cross-species comparison and consideration of probe performance, we selected array probes that can be reliably used for the crab-eating macaque genome. A DNA methylation assay using an EPIC array was performed on genomic DNA extracted from the brains of five crab-eating macaques. The obtained DNA methylation data were compared with a publicly available dataset. RESULTS: Among the 865 918 probes in the EPIC array, a total of 183 509 probes (21.2%) were selected as high-confidence array probes in the crab-eating macaque. Subsequent comparisons revealed that the data from these probes showed good concordance with other DNA methylation datasets of the crab-eating macaque. CONCLUSION: The selected high-confidence array probes would be useful for high-throughput DNA methylation assays of the crab-eating macaque.

Promoter Activity-Based Case-Control Association Study on SLC6A4 Highlighting Hypermethylation and Altered Amygdala Volume in Male Patients With Schizophrenia.

Associations between altered DNA methylation of the serotonin transporter (5-HTT)-encoding gene SLC6A4 and early life adversity, mood and anxiety disorders, and amygdala reactivity have been reported. However, few studies have examined epigenetic alterations of SLC6A4 in schizophrenia (SZ). We examined CpG sites of SLC6A4, whose DNA methylation levels have been reported to be altered in bipolar disorder, using 3 independent cohorts of patients with SZ and age-matched controls. We found significant hypermethylation of a CpG site in SLC6A4 in male patients with SZ in all 3 cohorts. We showed that chronic administration of risperidone did not affect the DNA methylation status at this CpG site using common marmosets, and that in vitro DNA methylation at this CpG site diminished the promoter activity of SLC6A4. We then genotyped the 5-HTT-linked polymorphic region (5-HTTLPR) and investigated the relationship among 5-HTTLPR, DNA methylation, and amygdala volume using brain imaging data. We found that patients harboring low-activity 5-HTTLPR alleles showed hypermethylation and they showed a negative correlation between DNA methylation levels and left amygdala volumes. These results suggest that hypermethylation of the CpG site in SLC6A4 is involved in the pathophysiology of SZ, especially in male patients harboring low-activity 5-HTTLPR alleles.

Global DNA hypomethylation and its correlation to the betaine level in peripheral blood of patients with schizophrenia.

Accumulating evidence suggests that aberrant epigenetic regulation is involved in the pathophysiology of major psychiatric disorders such as schizophrenia (SZ) and bipolar disorder (BD). We previously showed that the plasma level of betaine (N,N,N-trimethylglycine), a methyl-group donor, was significantly decreased in patients with first episode schizophrenia (FESZ). In this study, we identified decrease of global DNA methylation level in FESZ (N = 24 patients vs N = 42 controls), and found that global DNA methylation level was inversely correlated with scores on the global assessment of functioning (GAF) scale, and positively correlated with plasma betaine level. Notably, correlations between levels of betaine and its metabolites (N,N-dimethylglycine and sarcosine, N-methylglycine) were lower or lost in FESZ plasma, but remained high in controls. We further examined global DNA methylation levels in patients with chronic SZ (N = 388) and BD (N = 414) as well as controls (N = 430), and confirmed significant hypomethylation and decreased betaine level in SZ. We also found that patients with BD type I, but not those with BD type II, showed significant global hypomethylation. These results suggest that global hypomethylation associated with decreased betaine level in blood cells is common to SZ and BD, and may reflect common pathophysiology such as psychotic symptoms.

Altered expression of microRNA-223 in the plasma of patients with first-episode schizophrenia and its possible relation to neuronal migration-related genes.

Recent studies have shown that microRNAs (miRNAs) play a role as regulators of neurodevelopment by modulating gene expression. Altered miRNA expression has been reported in various psychiatric disorders, including schizophrenia. However, the changes in the miRNA expression profile that occur during the initial stage of schizophrenia have not been fully investigated. To explore the global alterations in miRNA expression profiles that may be associated with the onset of schizophrenia, we first profiled miRNA expression in plasma from 17 patients with first-episode schizophrenia and 17 healthy controls using microarray analysis. Among the miRNAs that showed robust changes, the elevated expression of has-miR-223-3p (miR-223) was validated via quantitative reverse transcription-polymerase chain reaction (qRT-PCR) using another independent sample set of 21 schizophrenia patients and 21 controls. To identify the putative targets of miR-223, we conducted a genome-wide gene expression analysis in neuronally differentiated SK-N-SH cells with stable miR-223 overexpression and an in silico analysis. We found that the mRNA expression levels of four genes related to the cytoskeleton or cell migration were significantly downregulated in miR-223-overexpressing cells, possibly due to interactions with miR-223. The in silico analysis suggested the presence of miR-223 target sites in these four genes. Lastly, a luciferase assay confirmed that miR-223 directly interacted with the 3' untranslated regions (UTRs) of all four genes. Our results reveal an increase in miR-223 in plasma during both the first episode and the later stage of schizophrenia, which may affect the expression of cell migration-related genes targeted by miR-223.

Evaluation of the usefulness of saliva for DNA methylation analysis in cohort studies.

INTRODUCTION: Epigenetic information such as DNA methylation is a useful biomarker that reflects complex gene-environmental interaction. Peripheral tissues such as blood and saliva are commonly collected as the source of genomic DNA in cohort studies. Epigenetic studies mainly use blood, while a few studies have addressed the epigenetic characteristics of saliva. METHODS: The effects of methods for DNA extraction and purification from saliva on DNA methylation were surveyed using Illumina Infinium HumanMethylation450 BeadChip. Using 386 661 probes, DNA methylation differences between blood and saliva from 22 healthy volunteers, and their functional and structural characteristics were examined. CpG sites with DNA methylation levels showing large interindividual variations in blood were evaluated using saliva DNA methylation profiles. RESULTS: Genomic DNA prepared by simplified protocol from saliva showed a similar quality DNA methylation profile to that derived from the manufacturer provided protocol. Consistent with previous studies, the DNA methylation profiles of blood and saliva showed high correlations. Blood showed 1,514 hypomethylated and 2099 hypermethylated probes, suggesting source-dependent DNA methylation patterns. CpG sites with large methylation difference between the two sources were underrepresented in the promoter regions and enriched within gene bodies. CpG sites with large interindividual methylation variations in blood also showed considerable variations in saliva. CONCLUSION: In addition to high correlation in DNA methylation profiles, CpG sites showing large interindividual DNA methylation differences were similar between blood and saliva, ensuring saliva could be a suitable alternative source for genomic DNA in cohort studies. Consideration of source-dependent DNA methylation differences will, however, be necessary.

Population-neuroscience study of the Tokyo TEEN Cohort (pn-TTC): Cohort longitudinal study to explore the neurobiological substrates of adolescent psychological and behavioral development.

AIM: Adolescence is a crucial stage of psychological development and is critically vulnerable to the onset of psychopathology. Our understanding of how the maturation of endocrine, epigenetics, and brain circuit may underlie psychological development in adolescence, however, has not been integrated. Here, we introduce our research project, the population-neuroscience study of the Tokyo TEEN Cohort (pn-TTC), a longitudinal study to explore the neurobiological substrates of development during adolescence. METHODS: Participants in the first wave of the pn-TTC (pn-TTC-1) study were recruited from those of the TTC study, a large-scale epidemiological survey in which 3171 parent-adolescent pairs were recruited from the general population. Participants underwent psychological, cognitive, sociological, and physical assessment. Moreover, adolescents and their parents underwent magnetic resonance imaging (MRI; structural MRI, resting-state functional MRI, and magnetic resonance spectroscopy), and adolescents provided saliva samples for hormone analysis and for DNA analysis including epigenetics. Furthermore, the second wave (pn-TTC-2) followed similar methods as in the first wave. RESULTS: A total of 301 parent-adolescent pairs participated in the pn-TTC-1 study. Moreover, 281 adolescents participated in the pn-TTC-2 study, 238 of whom were recruited from the pn-TTC-1 sample. The instruction for data request is available at: http://value.umin.jp/data-resource.html. CONCLUSION: The pn-TTC project is a large-scale and population-neuroscience-based survey with a plan of longitudinal biennial follow up. Through this approach we seek to elucidate adolescent developmental mechanisms according to biopsychosocial models. This current biomarker research project, using minimally biased samples recruited from the general population, has the potential to expand the new research field of population neuroscience.

Somatic mutations in the human brain: implications for psychiatric research.

Psychiatric disorders such as schizophrenia and bipolar disorder are caused by complex gene-environment interactions. While recent advances in genomic technologies have enabled the identification of several risk variants for psychiatric conditions, including single-nucleotide variants and copy-number variations, these factors can explain only a portion of the liability to these disorders. Although non-inherited factors had previously been attributed to environmental causes, recent genomic analyses have demonstrated that de novo mutations are among the main non-inherited risk factors for several psychiatric conditions. Somatic mutations in the brain may also explain how stochastic developmental events and environmental insults confer risk for a psychiatric disorder following fertilization. Here, we review evidence regarding somatic mutations in the brains of individuals with and without neuropsychiatric diseases. We further discuss the potential biological mechanisms underlying somatic mutations in the brain as well as the technical issues associated with the detection of somatic mutations in psychiatric research.